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Abstract 

During the last decades, the formalism of computational complexity has 
played a central role in advancing the understanding of the physics of spin 
glasses and in explaining the difficulty of simulating problems in statistical 
physics. The advent of quantum complexity theory has put us into the unique 
position to foster a similar understanding for systems of interacting quantum 
particles. The central problem in quantum mechanics, and more specifically in 
quantum chemistry and solid state physics, is to solve the Schrodinger equation 
for a system of electrons interacting via the Coulomb potential in an external 
field. We show that this problem is QMA (Quantum Merlin Arthur) hard by 
reducing it to a 2D lattice of qubits with a fixed Heisenberg interaction and 
local magnetic fields, which in turn can be reduced to a known QMA-complete 
problem. On our way, we find that the Hubbard model in the presence of 
local magnetic fields, one of the most intensively studied lattice models for 
fermions, is QMA-complete. 

This result is highly relevant within the context of Density Functional The- 
ory (DFT) , the most widely used and successful method for simulating systems 
of interacting electrons. In a series of papers that seem to be the most cited 
ones ever in physics, Kohn and collaborators proved the existence of a problem- 
independent black box universal functional [5] whose usage would considerably 
simplify the simulation of interacting electrons. The central problem in DFT 
has been to find better and better approximations to this universal functional. 
We show that the existence of good approximations would place the problem 
of interacting electrons in NP, implying NP=QMA. This provides strong ev- 
idence against the existence of an efficiently computable universal functional 
and poses bounds on its approximability. We finally discuss another widely 
used simulation method for interacting electrons, the Hartree-Fock ansatz, and 
show that finding the optimal energy there is NP-complete. 



1 Motivation and problem 



In this paper, we consider the problem of finding the ground state of a system of inter- 
acting electrons moving in an external potential from the point of view of computational 
complexity. This problem is encountered virtually everywhere in quantum chemistry as 
well as in condensed matter physics. For instance, the spatial configuration of a molecule 
is the one for which the energy of the interacting electrons moving in the nuclear po- 
tential, together with the electrostatic energy of the nuclei, gets minimal. Similarly, a 
rich variety of phenomena in solid state physics, in particular conductance and magnetic 
phenomena, can be understood by considering electrons moving in the periodic lattice 
potential, including such exciting phenomena as high-temperature superconductivity 
and the fractional quantum Hall effect. 

Let us first define the problem under consideration. We are given a system of iV 
electrons, which is described by a complex-valued £2 function ^(x\, . . . ,xn), where each 
Xi = (rj, Si) encodes the spatial coordinate ri € M 3 and the spin Sj G {0, 1}, and where 
the fermionic nature of the electrons requires ^ to be antisymmetric, 

*(..., Xi,...,Xj,...) = -V(...,Xj,...,Xi,...) 

for any pair of indices i and j. On the space of completely antisymmetric states, we 
want to find the ground state energy - the lowest eigenvalue - of the Hamiltonian 

ff = -^£ A * + £ur^H+£^)> (i) 
y * ' » « ' 

where the external potential V contains an electrostatic part (j>{ri) and a magnetic part 
B(?i), such that V(ri,S{) = 0(rj) — B(ri) ■ c?j with the Pauli matrices 5; L = (of , of , of) 
acting on the spin degree of freedom, and where the Laplace operator Aj is with respect 
to the i'th spatial coordinate. 

The computational task we consider is finding the ground state energy of this system, 
up to a polynomial precision in energy. More precisely, we will focus on the correspond- 
ing decision problem interacting electrons, where we have to decide whether the 
ground state energy is below some value a or above b, where b — a > l/poly(A r ) is the 
desired precision. 

The natural complexity class for ground state problems is QMA, Quantum Merlin 
Arthur [10]. A language L is in QMA if there exists a uniform family of polynomial-time 
quantum circuits V x such that for every x £ L, there is a polynomial-size quantum state 
\ip) (the proof) such that V x will accept \ip) with probability at least 2/3, while for 
x L, V x will accept any proof with probability at most 1/3. The canonical complete 
problem for the class QMA is LOCAL hamiltonian [10], where a local Hamiltonian on 
qubits is given and the problem is to decide whether its ground state energy is below a 
or above b for b — a > l/poly(A r ). Here, by local we mean interactions involving at most 
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a constant number of qubits each. Subsequently, there have been found more restricted 
ground state problems which are still QMA-complete, most notably for qubits on a 2D 
square lattice with a particular class of nearest neighbor couplings [13] and a ID chain 
of 12-level systems with nearest neighbor couplings [1]. 

In this article, we determine the computational complexity of INTERACTING ELEC- 
TRONS, which we show to be QMA-hard. Therefore, we consider the Hubbard model 
which describes electrons hopping on a lattice and is well-known to be a limiting case of 
the general problem of interaction electrons, Eq. (1); this reduces the original problem, 
defined on a continuous system, to a lattice problem for fermions. In a second step, we 
reduce the fermionic system to a system of distingushable spins, i.e. qubits, by showing 
that there is a domain in which the Hubbard model mimics the behaviour of qubits 
coupled via a fixed-strength Heisenberg interaction and subject to local magnetic fields. 
We then show how this Hamiltonian can be reduced to a model which is known to be 
QMA-complete, namely a two-dimensional square lattice of qubits coupled via so-called 
Pauli interactions [13]. 

Having shown QMA-completeness of interacting electrons, we turn our atten- 
tion to Density Functional Theory (DFT). DFT is the most commonly used method 
for finding the ground state energy of interacting electrons, and the most successful 
one. The central idea is to reduce the complicated many-particle wavefunction to the 
single-electron density, and encode all the interaction terms in a so-called "universal 
functional". This is motivated by the observation that the only variable term in the 
problem is the external potential which only depends on the single-electron density, 
while both the kinetic (i.e. phase-sensitive) and the interaction term are fixed, making 
the functional independent of the instance of the problem and thus "universal". The 
central task in DFT is now to find the universal functional, or a good approximation 
thereof. However, as we will show, the existence of an efficiently computable universal 
functional - or even of polynomial accuracy approximations - allows to place inter- 
acting electrons in NP, and would thus imply NP = QMA which is both considered 
unlikely and would resolve long-standing open problems, explaining the difficulty - or 
impossibility - of finding the universal functional. 

Finally, we turn towards Hartree-Fock, a precursor of DFT and still widely used 
in numerical calculations. Different from DFT, it is based on an ansatz for the wave 
function and is therefore not guaranteed to give the optimal ground state, while on the 
other hand it can be applied without having to guess some functional. We show that 
these differences also affect the computational complexity of the problem and prove that 
finding the minimum in Hartree-Fock is NP-complete. 

Our proofs will make excessive use of perturbation theory constructions. Rigorous 
derivations of perturbation expansions suitable for our purposes can be found in [9, 13], 
and we tacitly use these techniques in our perturbation expansions. 
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2 Reduction to a spin model 

The first part of proving the QMA-hardness of interacting electrons consists of 
reducing the complicated continuous three-dimensional many-electron system to a two- 
dimensional square lattice of qubits interacting via a fixed exchange interaction and 
subject to a tunable local field. 

We start by going to the two-dimensional Hubbard model [6], which describes a sys- 
tem of Fermions hopping on a lattice. It is a well-known special case of the Hamiltonian 
(1), obtained by putting a strong trapping field in the z direction and a periodic lattice 
potential in the x-y direction. This confines the electrons to localized fermionic modes 
on a two-dimensional lattice, described by fermionic creation operators a\ s . Here, i 
denotes the site on the 2D lattice, and s £ {0, 1} the spin state. 

The low-energy sector of the resulting system is described by the effective Hamilto- 
nian 



where <i,j> denotes nearest neighbors on the 2D lattice, n = a) a is the particle number 
operator, and <7j is the vector of Pauli matrices acting on the spin degree of freedom at 
site i. The first term in the Hamiltonian describes an electron tunneling from one site 
to the adjacent one without changing its spin; its strength is governed by the height and 
the width of the potential barrier. The second term gives the on-site repulsion: although 
each site only allows for one spatial mode, this mode can be occupied by one spin up and 
one spin down electron in a singlet state. However, in this case the energy is increased 
due to the Coulomb repulsion between the electrons, giving rise to the on-site repulsion 
U ; its strength depends on the width of the potential well. Finally, the rightmost term 
contains the contribution from the magnetic field which imposes a local field on the spin 
at each site i. Note that this is the only term which we need to tune locally. 

The Hubbard Hamiltonian is thus obtained from the perturbation expansion of the 
low-energy spectrum of interacting electrons in a particular external potential; a detailed 
derivation, including bounds on the higher-order terms, can be found e.g. in [2, 8]. Ex- 
cept for being a special case of interacting electrons, the Hubbard model is of high 
interest on its own, as it is believed to be the appropriate model for understanding ef- 
fects like high-temperature superconductivity, quantum magnetism, and heavy fermions. 
Indeed, it is one of the most intensively studied models in solid state physics, making 
the investigation of its computational complexity interesting on its own. A particular 
nice thing about the Hubbard model is that deciding its ground state energy is clearly 
inside QMA, since it can be mapped to a spin system via the Jordan- Wigner transform. 
Although this makes the Hamiltonian non-local, its ground state energy can still be 
measured efficiently [12], and with what follows, we have altogether that the Hubbard 
model is QMA-complete. 

To show the QMA-hardness of the problem, we transform the Hubbard Hamiltonian 
Eq. (2) to a Hamiltonian on a 2D square lattice of spins. Therefore, we choose the on- 





<i,j>,s 
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site repulsion energy U very large as compared to t (explicit choices for the scalings are 
given later), and operate the system in the so-called half-occupancy regime, i.e. where 
there are as many electrons as sites. Due to the large on-site repulsion, the hopping 
is suppressed to order t/U, and in the ground state each mode will be occupied by 
exactly one electron, providing a spin degree of freedom which we can use as a qubit. 
The coupling between the spins is obtained in second-order perturbation theory and 
is due to a virtual process where one spin tunnels to the adjacent site, interacts with 
the other spin, and tunnels back. Since two spins can only be at the same site if they 
form a singlet, this process decreases the energy of a singlet state by At 2 /U, and since 
(1 — ffiff j)/ 4 is the projector onto the singlet, the effective interaction Hamiltonian is 

H cxc = ^ ■ Bj - l - + 0(t 3 /U 2 ) . 
Omitting the energy offset, we are left with an effective lattice Hamiltonian 

#Hcis = \jY^B i -Bj-Y i § i -3 i + 0(Nt 3 /U 2 ) . (3) 

<i,j> i 

3 QMA-completeness of Hheis 

In the following, we show that deciding the ground state energy of the Hamiltonian 
i^Hcis with arbitrary local field is a QMA-complete problem. Note that it is clearly in 
QMA, as it is a special instance of the LOCAL HAMILTONIAN problem [10]. To prove 
QMA-hardness, we show that the problem can be reduced to the ground state problem 
of [13], where QMA-completeness is shown for a Hamiltonian on a 2D square lattice with 
interactions of the form 

<i,3> 

where and are Pauli matrices. Following [13], we call interactions of the form 
Ajj-A(jj) <g> Pauli interactions. 

We will now show how a 2D lattice with arbitrary Pauli interactions can be embedded 
into a 2D lattice with fixed Heisenberg interactions and local fields, Eq. (3). Therefore, 
we will use a number of so-called perturbation gadgets. These basic building blocks 
consist of a few qubits where a strong local field is applied to one of them, so that the 
interactions of the Hamiltonian induce an effective coupling on the remaining qubits in 
the low-energy subsector. An important observation made in [13] is that these gadgets 
can be combined: firstly, one can use several gadgets in parallel, as long as they do not 
share an interaction term, and secondly, one can stack gadgets atop of each other, as 
long as all new terms together are still in the perturbation regime, i.e. much smaller 
than the strong fields of the previous gadgets. In particular, this can be accomplished 
with polynomial-strength fields if the number of "gadget layers" is kept constant. For 
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Figure 1: Gadgets used to go from the 2D lattice with fixed Heisenberg inter- 
actions and local field to the 2D lattice with tunable Pauli interactions. 



the moment, we will discuss our gadgets one after another and postpone the discussion 
on how to appropriately connect them. 

We will use a total of 4 different gadgets which we apply atop of each other. The first 
gadget erases vertices from the 2D lattice, and we use it to create a sparse lattice where 
each effective edge consists of 8 edges with a Heisenberg coupling, see Fig. la. The second 
gadget reduces the full Heisenberg coupling to an XX coupling, thus mapping two edges 
to one (Fig. lb). The third gadget reduces two XX edges to one with an Ising-type 
coupling (Fig. lc), and the last gadget combines two Ising edges to obtain an arbitrary 
Pauli coupling of variable strength (Fig. Id). By properly concatenating the latter three, 
one can implement any Pauli coupling on each effective edge of the sparse lattice. The 
following derivations will make use of the techniques for perturbation gadgets as found 
in [9, 13]. 

The first gadget erases a node from the graph. Therefore, we apply a strong field B e 
to this qubit in, say, the z direction. The Hamiltonian is 

H CTasc = —B e al + \ °c-&q, 

q=l,r,t,b 

where in the following A = t 2 /U [cf. Eq. (3)], and where c, I, r, t, and b denote the 
center, left, right, top, and bottom qubit, respectively. Doing a perturbation expansion 
in \/B e , we find 




<8> |0)(0| c . 



Thus, the gadget induces an effective local Hamiltonian on the adjacent sites while 
decoupling them up to a perturbation of order X 2 /B e . Since the induced local terms 
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can be easily compensated by changing the corresponding local fields, the effect of this 
gadget is to erase the central qubit for sufficiently large B e . 

We now use the erasure gadget to cut a sublattice into the square lattice. Therefore, 
we apply the field B e to all sites we want to remove, so that we are left with a sparse 
lattice as in Fig. la, where each edge contains 8 edges of the original lattice. These edges 
will be used to build a single edge with an arbitrary Pauli coupling, using a sequence 
of gadgets each of which maps two edges to one while creating a new interaction. The 
first of these gadgets reduces the Heisenberg interaction to an XX interaction, i.e. an 
interaction of the form erf a? + cfej. Therefore, consider three qubits on a line and 
apply a field Bxx to the central qubit, 

H X X = —BxX&c + ^ ^ &c • &q ■ 
q=l,r 

Doing the perturbation expansion to second order, we obtain 

\2 / \3 



Hxx 



loc 



4B XX 



(of of + ofo y r ) + O 



B 2 
n xx 



|0)(0| c 



where H\ oc = — Bxx + A Yl q =i r i^q + (X/4Bxx)(Z q — 1)) collects all local terms, which 
can be compensated easily. Clearly, by choosing a different direction for Bxx, we can 
equally get couplings in x-z or in y-z. 

The next gadget reduces the XX coupling to an Ising type coupling crfaf- Therefore, 
we take a chain of three qubits coupled by the XX interaction (w.l.o.g. in x and z 
direction), and apply a field in Z direction to the second qubit, 

Rising = -Biol - n ( a l a r + <rf<4) 

q=l,r 

with /j, = X 2 /4:Bxx, which perturbatively gives 



sing 



-Bj 



Bj 



(of a? + 1) + O 



|0)(0| c 



and is thus locally equivalent to an Ising coupling. 

The last gadget is being used to get an arbitrary Pauli interaction with tunable 
strength. In order to get e.g. a coupling of the form of of, we take a three-qubit chain 
with two different Ising couplings and apply a strong diagonal field on the central qubit, 

tfpauli = -BpO^ - d{ofof + Olof) , 



with d = /U 2 /B/, and where o xy = \+) xy {+\xy ~ \-)xy{-\xy with \±) xy = (|0) + 
Vi\l))/V2, \/i = (1 + i)/y/2 vectors along the x + y direction. The resulting inter- 
action is is obtained from the perturbation series as 

-Bp - ^={of + of) - £-(of*v + 1 



Pauli 



2B P 



(4) 
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thus enabling us to construct any desired Pauli coupling. Note that this construction 
also allows to keep an Ising interaction af a* by applying Bp in x direction. 

The strength of the Pauli couplings can be tuned by appropriately changing the 
corresponding field Bp. However, there are a two caveats: On the one side, we want 
the 0(?9 3 / Bp) term above to vanish polynomially (in particular faster than the desired 
accuracy in the ground state energy) and are thus restricted to interaction strengths 
d/Bp < l/poly(iV). However, since the interaction strength in the QMA-complete 
construction for 2D lattices with Pauli interactions in [13] is polynomially bounded, we 
can rescale all interactions as well as the desired energy precision by a polynomial. 

On the other side, in order to obtain very weak couplings we would need to choose Bp 
arbitrarily large. Since, however, we want to keep all our Hamiltonian terms bounded 
by some polynomial in the systems size N, we cannot obtain couplings smaller than 
l/poly(iV). Again, this can be resolved by setting these coupling terms to zero, which 
results in an overall error of at most N/po\y(N) in energy. By choosing this sufficiently 
smaller than the desired precision, we can also avoid problems in this case. Clearly, 
setting the interaction terms to zero can be done by erasing those edges right in the first 
step when the erasure gadget is applied. 

Let us now go through the whole procedure of rescaling the interactions, choosing 
the local field strength, and assembling the gadgets, as to see that we stay in a poly- 
nomial regime. Recall that that when assembling the gadgets, the energy scale of all 
newly added gadgets should still be small compared to the large field underlying the 
perturbation expansion in the previous stage. We start from a 2D lattice with arbitrary 
Pauli interactions, where the interaction strength is bounded by a polynomial a(N), and 
with a desired energy accuracy l/b(N) with b a polynomial. Now choose a polynomial 
c 3> a 2 b and rescale everything by 1/c, thus obtaining a maximum interaction strength 
a/c and a precision 1/bc. 1 By omitting all interactions weaker than 0(l/N 2 b(N)c(N)), 
we make an error of 0(l/Nb(N)c(N)) in energy which is vanishing relative to the pre- 
cision. 

We now let U = 1 and t = 0(1/ Vbc) in the Hubbard model (2), leading to A = 
0(l/bc). Then, we choose the following magnetic fields: 

B e = O (Va(6c) 2 / 7 ) , 
Bxx = O (iV 4 A(6 C ) 2 / 7 ) , 
Bj = O (iV 2 A(6c) 2 / 7 ) , 
Bp = O (A(6c) 2 / 7 ) , 

which ensures four things: firstly, each field is much stronger than all following contri- 
butions, second, each field is much stronger than the corresponding perturbation A, fi, 

1 For convenience, we omit the parameter N when appropriate. 
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or thirdly, the permissible interaction strength $ 2 /2Bp is large enough to allow for 
0(l/N 2 b(N)c(N)) couplings, and finally, all the error terms from the perturbation ex- 
pansions are 0(N/ (fee) 9 / 7 ) and thus for sufficiently large c much smaller than the desired 
precision which is 0(1 /be). 

Thus, by assembling all the gadgets together we finally find that INTERACTING ELEC- 
TRONS is QMA-hard. 

4 Density Functional Theory 

As we have seen, finding the ground state energy of systems of interacting electrons is 
both central in physics and a computationally hard problem. Therefore, approximation 
methods have to be used in order to obtain results for realistic scenarios. One of the 
most popular methods for this problem is Density Functional Theory (DFT) [5, 14]. 
It is based on the observation that the only varying part of the Hamiltonian (1) is 
the potential term V which only depends on the single electron density, while all the 
other terms, and in particular the interactions, are fixed. Therefore, one can split the 
minimization over the set of antisymmetric A^-particle wavefunctions ^ in two parts: 
One considers the one-particle reduced density matrix given by 

P( r ) = / dr 2 ---dr N i $(r 1 ,s 1 ,...,r N ,s N )x 

, JM.3N-3 
s 1 ,s 1 ,s 2 ,---,s N 

^*(r 1 ,s' 1 ,r 2 ,S2, ■ ■ . ,r A r,s A r)|si)(s / 1 | 

which already fully determines the energy contribution coming from the the potential 
V. (Note that we can replace the spin density matrix by a vector in the Bloch sphere, 
making the symmetry between spatial and spin density more clear.) On the other hand, 
for a given p(r) one can minimize the energy of the fixed terms T + I (the kinetic and the 
interaction term) over all A^-electron states compatible with the marginal p[r) without 
having knowledge about V. Thus, the ground state energy of the problem is given by 

E = mm{tv(Vp) + F N [p]} , (5) 

p 

where the "universal functional" 

F N [p] = min (^Jy*(T + 

minimizes the fixed kinetic and interaction term over all A^-electron states compatible 
with p. Therefore, the central problem in DFT is to obtain the functional F/v, or at 
least a good enough approximation for the setup under investigation. 

In the following, we show that one cannot hope to obtain an efficiently computable 
universal functional F/v. Therefore, consider again the 2D Hubbard model which we 
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know to be QMA-complete, and let us show that we can specify p{r) for the ground 
state classically. 2 Therefore, note that we can specify the probability distribution over 
the discrete set of N lattice sites efficiently, and since the wavefunctions corresponding 
to the localized modes aj are known, this suffices to reconstruct the single-electron 
density matrix p{r). Thus, if Fn could be computed efficiently, one could provide an 
efficiently checkable witness for the ground state, and thus the ground state energy of the 
Hubbard model could be decided in NP, implying NP=QMA. Therefore, there cannot 
be an efficiently computable universal functional unless NP=QMA, which is considered 
unlikely. However, even if one is willing to accept NP=QMA, our result still emphasizes 
the hardness of finding the functional, since NP = QMA =>- NP = MA NEXP / 
MA NEXP <f_ P/poly [7], and whether NEXP contains hard Boolean functions is a 
long-standing open problem in computational complexity. 

Let us note that the argument for the impossibility of finding a universal functional 
also holds for approximations up to polynomial precision: for everything we have to 
know in the definition of QMA is the ground state energy up to a polynomial accuracy. 

One of the precursors of DFT which is still widely used is the Hartree-Fock method 
[11]. It is similar to DFT in that it reduces the iV-electron equation to a problem 
of individual electrons moving in an external field which depends only on the average 
electron distribution. Different from DFT, Hartree-Fock is based on a particular ansatz 
for the wave function and is therefore not guaranteed to give the true ground state 
energy; on the other hand, it is an iterative method which can be applied without prior 
assumptions, whereas in DFT some a priori guess on the form of the universal functional 
has to be made. 

The starting point is a Hamiltonian 

n = E H ij a hj + E H< S\i a H a ^ + • • • . ( 6 ) 

which is local, meaning the series terminates at some fixed H^ r \ and with a number 
of fermionic modes M > N; if it is derived from two-particle interactions as in the 
Schrodinger equation (1), r = 2. The Hartree-Fock method approximates the ground 
state of this Hamiltonian using the ansatz &]y---6||fi) with bi = J2 u ij a j (where 
is the vacuum). Note that this corresponds to an antisymmetrized product of single- 
particle wave functions, which is how Hartree-Fock is usually presented. 

In the following, we show that approximating the ground state energy using the 
Hartree-Fock method is an NP-complete problem. More precisely, we consider the prob- 
lem of deciding whether the lowest energy of (6) within the Hartree-Fock ansatz is below 
some a or above some b > a. We show that the problem is inside NP for up to an ex- 
ponential accuracy b — a and for any r, and that for NP-completeness a polynomial 
accuracy b — a < l/poly(iV) and r = 2 are sufficient. 

2 Note that this does not place the problem in QCMA [10] since we do not know how to obtain the 
optimal ^ from p even quantumly. 
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To see that the problem is in NP, note that a Hartree-Fock state is fully characterized 
by the u^'s, and that from there its energy can be computed efficiently. Conversely, the 
problem is shown to be NP-hard by mapping it to the ground state problem for Ising spin 
glasses which is known to be N P-hard [3] : Given an L x L x 2 lattice of two- level spins Si = 
±1 with a nearest neighbor Ising coupling TL = Yl JijSiSj, Jij G {0, —1, 1}, determine 
whether the ground state energy is the minimum one allowed by the individual Jjj's or 
not. Therefore, embed the TV = 11? classical spins into a fermionic system with IN 
modes occupied by N fermions. The modes come in pairs (an, o^i+i), and a Hamiltonian 
term Aw2i n 2i+i) A = 0{N 2 ) penalizes double occupancy, so that in the ground state 
exactly one mode per pair is occupied, giving an effective spin degree of freedom [12]; 
the coupling JijSiSj of these spins is realized as J2 pq =o x (— l) p+<? n2i+p"-2j+<?- As the 
ground state of the system is a classical spin state, it can be expressed as a Hartree-Fock 
state where bi = a<n or b{ = a2%+i, respectively, and since the classical Hamiltonian has 
a constant gap while perturbations from the penalized subspace are at most 0(1/A 2 ), a 
polynomial accuracy is sufficient to make the problem NP-hard. 

5 Outlook 

We have shown that interacting electrons, one of the core problems in physics and 
theoretical chemistry, is QMA-hard. On our way, we have proven QMA-completeness 
both for the Hubbard model and for qubits in a 2D lattice with a fixed Heisenberg 
coupling and local fields. We have then applied our findings to Density Functional 
Theory to show that the existence of an efficiently computable universal functional 
would place the Hubbard model and thus QMA in NP, making the existence of an 
explicit functional unlikely. Finally, we have shown that solving a system of interacting 
electrons using Hartree-Fock is NP-complete. 

An open question is whether interacting electrons can be placed inside QMA. 
The difficulty comes from the continuous nature of of the wavefunction, as compared to 
the Hubbard model where one can provide a proof in second quantization. Therefore, 
one would have to discretize the wavefunction in order to be able to encode the proof in 
a finite-size quantum state. In principle this should be feasible, since the kinetic term 
penalizes fast fluctuations in the wave function, thus justifying a coarse graining of the 
ground state witness. However, care has to be taken concerning the magnetic field since 
fast changes of the spin are not detected by the kinetic term. Therefore, in order to place 
interacting electrons in QMA, one will have to put constraints on the magnetic 
field bounding the speed of possible fluctuations. 

Another open question concerns the computational complexity of the problem of 
l/-representability [14]: given a local density p(r), does there exist a potential V(r) such 
that the Hamiltonian T + / + V [Eq. (1)] has a unique ground state with local density 
p(r)l The lattice version of the problem is the following: Given a local density p on 
the lattice, decide whether there exists a potential V such that H = Hq + V (for a 
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fixed Hq) has a unique ground state with density p. On the one hand, it can be proven 
in QMA whether a given density is the ground state to some local (spin or fermionic) 
Hamiltonian H - Merlin provides H, the ground state energy Eq, and an associated 
(possibly mixed) ground state 7, and Arthur either checks that 7 is the ground state of 
the local Hamiltonian H' = (H — Eq) 2 , or that it indeed gives rise to the proper reduced 
density. However, it is unclear how to simultaneously prove uniqueness of the ground 
state - it rather seems that disproving the uniqueness of the ground state for a given 
Hamiltonian is in QMA (Arthur provides two non-identical ground states). Note that 
without the requirement of uniqueness of the ground state - i.e. asking whether there 
exists a potential V such that some (mixed) ground state 7 gives rise to p - the lattice 
version of F-representability can probably be decided efficiently [4]. 
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